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Abstract 

New volume and volume-surface integral equations are presented 
for modeling inhomogeneous dielectric regions. The presented inte- 
gral equations result in more efficient numerical implementations and 
should therefore be useful in a variety of electromagnetic applications. 
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1 Introduction 


The modeling of inhomogeneous dielectrics via an integral equation ap- 
proach is traditionally accomplished via the introduction of equivalent vol- 
ume electric and magnetic currents [1] - [8]. For a dielectric with non-trivial 
permittivity and permeability this type of modeling implies six scalar un- 
knowns at each volume location. As a result, the implementation of the 
resulting integral equation is computationally intensive and has excessive 
storage requirements. 

In this paper it is demonstrated that any inhomogeneous dielectric mate- 
rial, regardless of its permittivity and permeability profile, can be modeled 
by a single electric or magnetic current density. Alternatively, either the 
electric or magnetic fields within the dielectric can be used as the unknown 
quantities. It appears though that one must pay a price for resorting to 
these reduced-unknown and/or kernal-singularity representations. Specifi- 
cally, because they involve derivatives of the unknown quantities, a higher 
(at least linear) basis function is required for discretizing the resulting inte- 
gral equations. However, it is possible to relax this requirement by resorting 
to a new volume- surface field representation. In this case, either the undif- 
ferentiated electric or magnetic field within the dielectric is the unknown 
quantity along with the corresponding tangential electric or magnetic fields 
on the outer boundary. Provided the dielectric volume is not composed of 
a single thin layer, this volume- surface integral equation still represents a 
nearly fifty percent reduction in the number of unknowns relative to tradi- 
tional implementations. 

2 Volume Representations 

Let us consider the dielectric/ferrite volume Vd, shown in Fig. 1, having 
relative constitutive parameters e r and fi r which are arbitrary functions of 
position. Assuming some exterior excitation, (E*,H‘), the total field can 
be written as 

E = E' + E s H = H‘ + H S (1) 

where (E s , H 5 ) are the scattered fields caused by the presence of the dielec- 
tric. Traditionally [1] the scattered fields are formulated in terms of the 
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equivalent currents 


Je, = jk 0 Y 0 (e r - 1)E, M eg = jk 0 Z 0 (fl T - 1)H (2) 

with k 0 and Z 0 = 1 /Y 0 being the free space wavenumber and intrinsic 
impedance, respectively. In terms of these effective or equivalent current 
densities, the scattered field is given by 

E s = JJJ v [V x r o (r, r') • M eg (r') + jk 0 Z o T o (r, r') • J e ,(r')] dv' (3) 

in which r and r' denote the observation and integration points, respec- 

tively, 

r o (r,r') = -fl+^lG.(r,r'), (4) 

is the free space dyadic Green’s function, 

V x F 0 (r, r') = -VGoir^xI, (5) 


Go( r,r') 


e -jfc 0 |r-r'| 
47r|r — r'| 


(6) 


I = xx + yy + zz 


is the unit dyad and H s is given by the dual of (3). By substituting (3) and 
its dual into (1) and then into (2), we obtain the coupled set of integral 
equations 


•Mr) 

jk 0 Y o (^t T 1 ) 


r€V d 


Meg(r) 

jWfir ~ 1 ) 


H" = H ! r € V d 


for a solution of the equivalent currents J e? and M e? . 

The aforementioned formulation appears to be the only approach that 
has so far been utilized for three dimensional implementations. However, 
as noted in the literature [4, 5, 9], the singularity of the kernal in (3) 
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presents numerical difficulties. Also, for non-trivial permeability, six scalar 
unknowns axe involved in the solution of (7). The first of these difficulties 
can be alleviated by resorting to higher order basis functions and expressing, 
for example, E 4 as 


E s 



x VG 0 (r, r') — j k 0 Z 0 3 eq G 0 (r, r') 


— • J e? (r , )VG 0 (r, r') 


dv' 


( 8 ) 


which is the volume equivalent to the Stratton-Chu surface integral equa- 
tion. Likewise, the scattered magnetic field to be substituted in (7b) can 
be expressed by the dual of (8). 

Although this appears to be the most popular approach in modeling 
three-dimensional dielectrics, it can be shown that there are several other 
ways to formulate the problem. Most importantly, it can also be shown 
that (7) can be replaced with an equivalent system which involves only 
three (not six) scalar unknowns. Specifically, from Maxwell’s equations [10] 
the radiation of M e , is indistinguishable from the radiation of the electric 
current 


j: 


eq 


V X Me, 
j k 0 Z 0 


( 9 ) 


This can be combined with (2) giving a single equivalent electric current 


J" = jkoYJer ~ 1)E + V x [(/x, - 1)H] 


( £ r - 1) 


V x H + V x [(/i r - 1)H] 


( 10 ) 


for representing the scattered fields (E 4 , H s ). From the dual of (3) we then 
obtain that the scattered magnetic field due to the current density (10) is 


H 5 


= /// v< K(r,r')xI] 


eJ $T v * H < r '> 


+ V' x [Mr') - l)H(r')] U' 


( 11 ) 
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in which V' implies differentiation with respect to the primed/integration 
coordinates. When this is used in (1) we deduce the integral equation 


H*(r) = H(r) - /y / [VG 0 ( r,r')xl 
+ V' x [(/z r (r') - l)H(r')] }du' 


l er(r') 


V' x H(r') 


reV d (12) 


where the unknown quantity is now the magnetic field within V d . Using a 
similar procedure it can be also shown that the scattered field may instead 
be represented by the radiation of a single magnetic current density 

M" = - (//r ~ 1) V x E - V x [(e r - 1)E] (13) 

fir 


From the first of (1) and (3), we then deduce the integral equation 



^r(r') - 1 

Mr') 


V' x E(r') 


+ V'x[(e r (r')-l)E(r')] W 


')]}' 


(14) 


which as expected is the dual of (12). We observe that the kernel singular- 
ity associated with (12) and (14) is the same as that associated with (8). 
In addition, as in the case of the integral equation (7) in conjunction with 
(8), linear expansion functions such as those in [3] or [4] are required for 
the discretizing (12) and (14). Thus, even though the new integral equa- 
tions (12) ,and (14) have half the unknowns, this was not achieved at the 
expense of increasing the kernel’s singularity or the order of the expansion 
basis required in their implementation. It is remarked that special forms 
of these integral equations have already been successfully implemented for 
two dimensional applications [11, 12]. 


3 Volume- Surface Representation 

The requirement to employ linear basis in connection with the implemen- 
tation of (12) and (14) can be relaxed by resorting to a volume- surface 
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integral equation (VSIE) such as that derived in [13] and [14] for two di- 
mensional simulations. To do so we begin with (3) which in conjunction 
with (2) can be rewritten as 

E- = E; + E‘„ = -kl J J Mr') - 1) E(r') • r„(r, r') 

-jk„Z,y xJJJ [MO - 1) H(r')G„(r, r')dv' (15) 

where is associated with the second integral and represents the field 
due to the magnetic equivalent current defined in (2). Setting 
H = V x E/ jk 0 Z 0 n r in this integral, and invoking the identities 

V x [V' x <f>E] = V x [VV x E] + V x [<j>V x E] 

V{<pxp) = ifiV <f> + fiV ip 

we obtain 

E;=F^+F" + F; 3 (16) 

with 

F ~ = Vx //l v ' x {( 1 -^) G ^ r ') E M} d *’' < 17 > 

F " - ( is > 

f ~ = v *JJL { G »< r - r '> v ' fcy * E(r ')} . < i9 > 

These integral expressions can be simplified through the use of various 
integral and differential identities. 

The volume integral in (17) can be transformed to a surface integral by 
invoking Stoke’s identity 

1 1 Iv X A ^ dV> = § s ^ h ' X A ^ l$ ' (20) 
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where Sd is the surface enclosing Vj and h' = n( r') denotes the outward 
unit normal to the surface Sd- We have 


= - §s. (l - [ft' x E(r')] x VG.(r, r')*' (21) 

which is an integral involving the undifferentiated tangential electric field 
over the surface enclosing Vd- Turning now to the integral in (18) we first 
rewrite it as 

K = - f [ [ - -Lrl V x [V'G„(r, r') x E(r')] dv’ (22) 

J J Jv d [ AM r ). 

and we note that [15, p. 487] 

V x [VG 0 x E(r')] = E(r')V 2 G 0 - E(r') • VVG 0 (23) 

Then, upon invoking the differential equation 

V 2 G 0 (r, r') + klG 0 ( r, r') = -6 ( r - r') (24) 


where <5(r') denotes the Dirac delta function, it follows that 

F " - JJL I 1 - ^?)) E < r ') • w + 1 1 - ^j] EM < 25 > 


Again, this involves only the undifferentiated electric field within the di- 
electric’s volume. Finally, the last integral in (16) can be readily simplified 
and written as 


F 


53 

m 


(26) 


( 



When (21), (25) and (26) are substituted into (16) and then into (15), 
we find that the total scattered field can be expressed as 

+//l VG ^ r ') x { v '[^] XE ( r ' ) } A ' 

- §s t [l - K X E(r')] x VG„(r,r')di' 

+ [ 1 "^] E(r) (27) 

For two dimensional simulations where the material parameters and the 
fields are invariant with respect to z, this expression can be readily shown 
to reduce to the VSIE given by Jin, etc. [13, equs. 28 and 31]. Expression 
(27) is also similar to the VSIE given by Tai [16]. However, Tai’s expression 
was left in terms of differentiated field quantities and is only applicable for 
homogeneous dielectrics. 

To obtain an integral equation on the basis of (27) we substitute this 
into the first of (1) and upon taking the principal value of the appropriate 
integrals we have 

+ jiLi E ^* v i^}} xvG ^ r 'w 

- §s,-s, [l - [ft' x E ( r ')l x VG °( r > r ') dv ' 

f E(r) r not in Vj 

+ E ‘ = | H 1 + sfo] E ( r ) r on s * ( 2S ) 

. 5 [ e r(r) + 2// r (r)] E(r) r in V d 


S 



In this, V 0 is a vanishingly small spherical volume whereas S 0 is a vanish- 
ingly small hemispherical surface both having their centers at r. As given, 
(28) can be discretized via the moment method or some other technique 
for a solution of E(r) within the dielectric. Its kernal has, of course, the 
same singularity as (7a) but involves only a single unknown vector field 
in comparison with the two vector unknowns appearing in (7). If linear 
rather than pulse basis are employed for the solution of (28), then it may 
be desirable to rewrite the first integral of (28) in the form given by (8) 
with M e , = 0 and 


Je</ 


jko e r (r)// r (r) 
Z 0 " f^r(r ) 



(29) 


However, in this case one could also resort to the alternative integral equa- 
tions (12) or (14). Of course, the dual of (28) is another integral equation. 
Further, linear combinations of (28) and its dual or (12) and (14) can be 
utilized if so desired. 

In closing, we remark that if // r and/or e T are discontinuous within Vj, 
the surface integral in (27) and its dual must then be replaced by 

H$s di [<( r ') - u-(r')] [n,(r') x F(r')] x VG 0 (r,r')ds' 

i 


where F = E or H. Here, Sj, denotes the ith discontinuous surface within 
Vd-, n,( r) is the unit normal to S, pointing from the — side to the -f side 
(outermost side) and n± denotes the inverse relative dielectric constant at 
the + or — side of the surface 5^, . In particular u' = 1 / fi' r for the E-field 
integral equation (27) and u' = 1/eJ. for the H-field integral equation. 


4 Conclusion 

Some alternative formulations were proposed for modeling three-dimensional 
inhomogeneous dielectrics. These are summarized in figure 2 and the aim 
of the investigation was to generate integral equations for the fields within 
the dielectric scatterer utilizing the minimum number of unknowns and the 
least singular kernels. A purely volume integral equation was derived in- 
volving half the unknowns required with traditional equations for ferrite 
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materials. The implementation of this reduced-unknown volume equation 
implies use of (at least) linear basis functions and to relax this requirement 
a volume-surface integral equation was derived. All of the integral equa- 
tions presented here appear to be more efficient than the traditional ones 
without compromising the kernel’s singularity. They should thus be found 
useful in a variety of radiation, scattering or SAR applications. 
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Fig. 1 . Illustration of the inhomogeneous dielectric volume V d 
enclosed by the surface S d . 



<c 


jk 0 Y 0 (e r -l)E jkoZo(M)H 






Fig. 2. 


Different volume equivalent currents for modeling the scattering 
by the inhomogeneous dielectric volume V d in figure 1 . 



